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ABSTRACT 

Ammonia stripping is a widely used process for removal of 
ammonical nitrogen from water. However, the expenditures 
involved, both in terms of capital, and operation and maintenance 
costs, are quite considerable. Thus even a few percent saving at 
the design stage would be of significant value. In the present 
study an attempt is made to develop a rational approach to 
minimum cost design of ammonia stripping towers. Procedures for 
optimal design of both counterflow and crossflow ammonia 
stripping towers are developed. The feasibility of preheating 
water and/or air for improving the efficiency of ammonia removal 
has also been investigated. Optimization problem is framed as 
minimisation of capital and operational cost subject to the 
constraints determined by process requirements. Interior penalty 
function method has been used as an optimization technique. 
Results show that in most of the cases crossflow towers are 
efficient than counterflow towers. The results also indicate that 
the option of preheating water and/or air to improve ammonia 
removal efficiency may be worth considering in some situations. 

KEYWORDS 

Gas transfer. Ammonia stripping. Packed bed towers. Counterflow 
towers, crossflow towers. Non - linear programming, optimization. 


Least cost design. 



INTRODUCTION 


Various compounds containing nitrogen are becoming increasingly 
important in wastewater management programs because of the many 
effects that nitrogenous materials can have on the environment. 
Nitrogen in its various forms, can deplete dissolved oxygen levels in 
receiving waters, stimulate aquatid growth, exhibit toxicity towards 
aquatic life, affect chlorine disinfection efficiency, present public 
health hazard, and affect suitability of water for reuse. Major 
processes for removing nitrogen from wastewater are biological 
nitrification-denitrification, break point chlorination, selective ion 
exchange for ammonia removal, and air stripping for ammonia removal 
(ammonia stripping). These are the processes which are technically and 
economically most viable at the present time. Out of these methods, 
ammonia stripping is one of the most economical means of nitrogen 
removal. The ammonia stripping concept is based on very simple 
principles. Due to its simplicity, it offers a reliable means of 
ammonia removal when applied under appropriate conditions. Ammonia 
stripping is carried out in stripping towers. There are two basic 
types of towers now being used in full scale operation: counter 
current and cross current towers. 

Design parameters for ammonia stripping towers include air to 
liquid ratio, tower depth and loading rates. Proper choice of these 
variables is required for efficient and economical nitrogen removal. 
Also the efficiency of ammonia stripping towers depend on the influent 
water and air temperature, with efficiency of ammonia removal reducing 
drastically at low temperatures. So an attempt should be made to study 
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the feasibility of preheating air and/or water for improved ammonia 
removal efficiency. 

Ammonia stripping towers are mostly designed using empirical 
design procedures. Performances of stripping towers have been 
documented from pilot plant studies. Pilot plant studies of ammonia 
stripping were conducted (Slechta and Culp, 1967 ; EPA, 1971 ) at South 
Tahoe public utility district. Based on these studies it was suggested 
that a liquid loading rate of 5000 kg/hr/m 2 is compatible with 
efficient tower operation at 6 to 7 m packing depth. Pilot studies 
were also conducted in Orange County water district (Uesner and Argo, 
1973) and they proposed an air to liquid ratio of 3000 m 3 /m 3 . Effect 
of temperature on ammonia removal was studied at the Blue Plains plant 
in Washington D.C. (O'Farrel et al . , 1973). It was reported that at 
10°C the maximum practical removal efficiency drops to about 75 
percent. Voluminous literature is also available on the practical 
aspects of the operation of the ammonia stripping towers (Culp and 
Culp, 1971; Gonzales and Culp, 1973; Kepple, 1973 and Culp, 1974 ). 
Based on various pilot plant studies Environmental Protection Agency, 
USA published a manual on nitrogen removal (EPA, 1975 ) which covers 
many aspects of the design of ammonia stripping towers. 

A review of the literature related to ammonia removal using 
packed bed towers reveals that many investigations have been carried 
out in the last two decades to improve upon the design procedures, 
practical aspects, and hence performance of packed bed towers. These 
investigations deal with mathematical modelling (Snow and Unek, 1969 ; 
Roesler et al . , 1970, 1971), development of computer programs for 
analysis of temperature profiles in the tower (Chiem et al . , 1981 ), 
improved mass transfer models (Rolf et al . , 1982 ), hydraulics of 
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packed bed towers (Billet, 1983), application of absorption data for 
design and scale up of packed bed towers (Billet and Hackowick, 1984), 
etc. An attempt on optimization of ammonia removal from wastewater was 
made by Rozhdov and Skidan (1980). A simulation model was described 
for optimization of removal of ammonia from wastewater by desorption 
of ammonia after alkali pretreatment, followed by oxidation of 
residual ammonia with active chlorine. This model makes it possible to 
determine the optimal values of such parameters as aqueous phase pH, 
air consumption, tower cross section, packed height of raschig ring 
and chlorine and lime dose. A model for optimization of packed tower 
based on optimal calculation of the gas flow pressure drop relative to 
a group of system variables was proposed by Rodriguez et al (1986). 
However, these investigations emphasised on maximizing ammonia removal 
efficiency instead of minimum cost design. The present study 
specifically focuses on this aspect of the ammonia striping towers. 

SCOPE 

The major objective of the present study is to evolve a procedure 
for minimum cost (optimal) design of ammonia stripping towers. An 
analysis of the design procedure (Roesler et al . , 1971) indicate that 

the equations involved are nonlinear in nature which necessitate 
application of nonlinear optimization techniques. Interior Penalty 
function method (Rao, 1978) which is a popular and versatile technique 
for nonlinear optimization is used. Specifically, the entire work was 
directed on the following lines: 

1. Development of analysis and design procedure for counterflow 
ammonia stripping towers. 

2. Development of analysis and design procedure for crossflow ammonia 
stripping towers. 
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3. Formulation of procedure for optimal design of counter and cross 
flow ammonia stripping towers with and without the option of 
preheating influent air and/or water. 

4. Evaluation of cost functions for various components influencing the 
process, including the cost of preheating water and/or air. 

5. Analysis of the optimization results. 

THEORETICAL BACKGROUND 

The process of ammonia stripping consists of transf ering ammonia 
gas dissolved in water to air which is in contact with water. The 
usefulness of ammonia stripping tower lies in the fact that it helps 
this transfer to occur in an efficient way. 

Equilibrium Relationship 

To study the process of ammonia stripping, it is convenient to 
first consider the equilibrium relationships. The relationship between 
concentration of ammonia in water and air can be depicted as shown in 
Fig. 1. The initial portion of this curve may be approximated by a 



Fig. 1. Equilibrium distribution of ammonia 
between a gas and a liquid phase 
at constant temperature . 
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straight line. If the slope of this line is at, then the concentration 
of ammonia in water and air can be related by the expression 


y Ai « a . x Ai (01) 

where, y^ is the equilibrium concentration of ammonia in air, x^. is 
the equilibrium concentration of ammonia in water, and a is a 
constant. This relationship is also known as Henry's law of diffusion. 
For dilute solutions (ammonia concentration < 5000 ag/1 , Treybal , 
1980) this law gives good approximation of the actual situation. Thus, 
in general, whenever a substance is dissolved between two insoluble 
phases, a dynamic equilibrium is established with equilibrium 


concentrations governed by Henry’s law. 

The curve shown in Fig. 1 is unique for a particular liquid-gas 
system and the substance to be transferd. Nevertheless the following 
principles are common to all systems involving distribution of a 
substance between two insoluble phases. 

1. For a fixed set of conditions, i.e. temperature and pressure, there 
exists an equilibrium relationship which may be shown graphically in 
the form of an equilibrium distribution curve by plotting the 
equilibrium concentrations in the two phases one against the other. 

2. For a system in equilibrium, there is no net diffusion of 
components between the two phases. 

3. For a system not in equilibrium, diffusion of components between 
the phases will occur in such a manner as to bring the system to a 
condition of equilibrium. If sufficient time is available, equilibrium 
concentrations will eventually prevail. 

Rate of Hass Transfer 

The mass transfer takes place according to the two film theory 
(Fig. 2). According to this theory there is no resistance to solute 
transfer across the interface separating the phases, and as a result 
the concentrations at the interface of the two phases are equilibrium 
concentrations. Mathematicaly , 
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Fig. 2. The two film theory concept . 

N. - F_ ln[ ---] * (02) 

1 " y AG 

where, is the rate of gas transfer/unit area, is the gas 

transfer coefficient, and y^ is the ammonia concentration in gaseous 
phase. Thus it is evident that for efficient mass transfer: 

1. The driving force i.e. the departure of the system from equilibrium 
should be more. 

2. The area provided for mass transfer should be more. 

The first condition is achieved by letting water and air mix in a 
certain way i.e. in counter or cross current mode. The second 
condition is achieved by packing the tower. Packing is usually 
provided to minimize the film resistance to gas transfer by continuous 
forming, splashing and reforming of drops. 

ANALYSIS OF COUNTERFLOU TOWERS 

The detailed analysis of counterflow ammonia stripping tower has 
been presented by Treybal (1980). According to this analysis, the 
height of the packed bed tower (Z) can be expressed as 
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H tG* N tG 


H 


tG 


Jr- 


initial 

dy 


(03) 


f Ai 


r final 


where, H tG is the height of transfer unit for ammonia transfer, N tG is 
the number of transfer units, and y initlal , y final and y are the 
ammonia concentrations in gas phase at the inlet, exit and at any 
level respectively in the tower. 

Solution Procedure 

The design of packed bed tower for ammonia removal requires 
evaluation of both temperature and ammonia concentration profiles in 
the tower. Distribution and rate of ammonia transfer between liquid 
and gaseous phase are dependent upon temperature and hence temperature 
profile in the tower has to be evaluated first. 

Computation of Temperature Profile 

The schematic diagram of the counter flow ammonia stripping tower 
is given in Fig. 3. The tower is divided into elements and each 
element is solved separately starting from bottom of the tower. Since 
the concentration of ammonia in both water and air is very less, it 
can be safely assumed that 


L 1 * L 2 


G , and 
L 


where, G^, and G are the total gas flow rates at the entrance, exit 
and at any level respectively in the tower, and L^, L 2 and L are the 
corresponding total liquid flow rates. The equations used for the 
computation of temperature profile are developed based on heat balance 
and heat transfer for an isolated element (Fig. 4), and can be written 
as follows (Me A dams , 1954). 



e 



x 10 y 10 


t i 

®1 l 2 

Fig. 3. Schematic representation of 
counterfiow ammonia stripping 
tower 
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<h j ' h j41 5 * L ’ c p CIWj - Tv J+1 ) COO 

Ch j ' h i*l ) ’ a/HT) <h ’cj ‘ h cj 5 A* 


where, G’ and L* are air and liquid mass flow rates respectively, HT 

is the height of transfer unit for cooling, h. and h , are the 

J c j 

enthalpies of air at the entrance and at the centre of the j th 

* 

element respectively, h . is the enthalpy of air at the centre of the 

^ J 

j th element corresponding to the temperature of water, C is the 

P 

specific heat of water. Tv ^ is the temperature of water at the 
entrance of the j th element, and / z is the incremental height. The 
value of HT can be determined empirically as mentioned by Roesler et 


al (1971). 

HT * (407.03/A p ) (L'/G’) 0,3 .....(06) 

Here, A p , the specific surface area of packing is given by 

A p - 60.96 (H t + V t )/(H s . Vg) (07) 


where, H^., V^, Hg and Vg are horizontal thickness of wood slats, 

vertical thickness of wood slats, horizontal spacing between wood 
slats and vertical spacing between wood slats used as packing 
respectively. Empirical regression relationships (Perry, 1963) are 


used to correlate wet-bulb temperature and enthalpy of air as 
h * 0.5542 [0. 00008 (1.8 Ta + 32) 3 - 0.0064 (1.8 Ta +32) 2 

+ 0.54789 (1.8 Ta + 32) - 1.62754] (08) 

Ta « 0 . 555 { [ 0.00016632 (h/0.5542) 3 - 0.036058 (h/0.5542) 2 

+ 3.1195 (h/0.5542) - 0.39736] - 32.0} (09) 


where, h is the enthalpy of air, and Ta is the temperature of air. 
The flow chart for computation of temperature profile is presented in 


Fig. 5. 


start 
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Fig. 5. Flow chart for computing temperature and ammonia concentration 
profiles in counterflow ammonia stripping tower. 
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Computation of Ammonia Concentration Profile 

1 

The equations used for the computation of ammonia concentration 
profile are developed based on mass balance and mass transfer 
relations with reference to Figs. 3 and 4. 

® ^ y j ~ ^ ~ *j+i^ (10) 

y j " y j + l * Cy c j “ y cj 5 L\ z C11) 

Eq. 11 is the difference form of Eq . 3. Here, y. and y . are ammonia 

J ^ J 

concentration in gaseous phase at the entrance and at the centre of 

the j th element respectively, y is the ammonia concentration in 

c j 

gaseous phase at the centre of j th element in equilibrium with 

ammonia concentration in water, and x. is the ammonia concentration in 

J 

liquid phase at the entrance of the j th element. The value of H^_ is 
given as follows (Roesler et al . , 1971). 

H tQ - 2.929 (G/L) 1 - 322 (12) 

The computational procedure for ammonia concentration profile is 
presented in Fig. 5. 

Computation of Air Pressure Drop 

The air pressure drop in a counterflow packed bed ammonia 
stripping tower is a function of tower height, air flow rate and the 
nature of packing. Generally an empirical correlation (Johnstone and 
Singh, 1937) is employed for the computation of pressre drop. 

/\P « 0.482 . f . z . (G’) 1,8 (13) 

where, /\ P is the air pressure drop in cm of water, and f is the 

fanning fraction factor (Table 1). 

ANALYSIS OF CROSSFLOW TOWERS 

The detailed analysis of crossflow ammonia stripping tower has 
been presented by Thibodeaux (1969). According to this analysis, the 
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Tabl* 1. Friction factors corresponding to various wood grid sisss. 

(H,p = 0.0635 eras ) 


H, 


f x 10 


8 


2.54 

5.08 

10.16 

20.32 

10.16 

20.32 

30.48 

10.16 

20.32 

30.48 

10.16 

20.32 

30.48 


1.5875 

1.5875 

1.5875 

1.5875 

3.1750 

3.1750 

3.1750 

4.4450 

4.4450 

4.4450 

5.7150 

5.7150 

5.7150 


18.90 

13.70 

8.30 

5.70 

2.43 

1.79 

1.20 

1.73 

1.16 

0.95 

0.87 

0.68 

0.63 


incremental width /_ w of a crossflow tower can be expressed as 



y i+l j 

where, is the ammonia concentration in gaseous phase at the 

entrance to the element (i j). 

Solution Procedure 

The schematic diagram of the crossflow tower is shown in Fig. 6. 
Similar to the analysis of counterflow tower, analysis of crossflow 
tower involve two steps, viz., calculation of temperature and ammonia 
concentration profiles. 

Computation of Temperature Profile 

The schematic diagram of the crossflow ammonia stripping tower is 
given in Fig. 7. The tower is divided into elements and each element 
is solved separately. The details of a typical element are shown in 
Fig. 8. The following two equations (He Adams , 1954) are used for the 
computation of the temperature profile. 




Fig. 6. Schematic representation of crossflow ammonia 
stripping tower . 



Fig. 7. Details of an element of crossflow ammonia 
stripping tower. 
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Fig. 8. Conceptual plan of crossflow ammonia 
stripping tower. 

g- ch itl j-n, p a* - <-• a-i J+1 - tvj j5 c ? a« (lg) 

h i j - h i+l j - < h \u - h cij> (“5 

where, h . . and h . . are the enthalpies of air at the entrance and at 

* J i J 

* 

the centre of an element (i j) respectively, h c . ^ is the enthalpy of 
air at the centre of element (i j) corresponding to the water 
temperature, and Tw^ is the temperature of water at the entrance of 
the element (i j). Eq. 15 is the difference form of Eq. 14. Substi- 
uting /\ z = /\w in Eq. 15, the following equation can be obtained 

©’ Ch i+1 j “ hi j) * L’ (T^i j+1 - Tw a 3 ) C p (1?) 

Empirical regression relationships given by Eqs. 8 and 9 are used to 
relate enthalpy and wet-bulb temperature of air. The algorithm for 
calculation of temperature profile is given in Fig. 9. 




Fig. 9. Flow chart for computing temperature and ammonia 

concentration profiles in crossflow ammonia stripping 
tower . 
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Computation of Ammonia Concentration Profile 

The equations used for computation of ammonia concentration 
profile are developed based on mass balance and mass transfer 
relations with reference to Figs. 7 and 8. 


(2.A/BB) + A - m 2 

*i j+i “ [ 3 *i i + I 3 i 

J (2.A/BB) + A + m 1 J (2.A/BB) + A + m L j 

where , 

m = 0.1117 e °* 02 * 12 ( 1 -8Tw i ^ + 32.0) 

A « L/G 

BB = ZiWH tG 

and 


(18) 

(19) 

( 20 ) 
( 21 ) 


V i j “ y i + l j “ L/G Cx i j + l - x i j> ( 22 3 

Here, x„ and y^ are the ammonia concentrations in liquid and gaseous 

phase at the entrance of element (i j) respectively. The algorithm for 

calculation of ammonia concentration profile is given in Fig. 9. 

Computation of Air Pressure Drop 

The air pressure drop in a crossflow tower is related to air flow 
rate, width of the tower, nature of packing and density of air. The 
relationships used in the present study are adopted from Kelly and 
Swenson (195€>). According to these relationships, the pressure drop 
consists of two components: 

1. the pressure drop caused by filling and support members. 

[j x * (4.77 . QL . G' 2 . U)/(H g P) (23) 

— 8 

where, (X is a constant equal to 0.26 x 10 , W is the width of an unit 

of a crossflow tower, and P is the density of air. 

2. pressure drop due to obstruction by water drops. 

I \? 2 = 2.24 x 10 -3 (.p. S F °* 5 . L*. G £ 2 . 2)/(Vg. p ) 


(24) 
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where, Q ia a constant equal to 0.07 x 10~ 8 , G is the equivalent air 
siass velocity for pressure drop relation given by 

G £ « Vt (0 - 0422 G ' 2 ) + (489.15 G*. p . S p °‘ 5 > + (1596029 .0 p 2 . S ) 

(25) 

and Sp is the vertical fall of a water drop used in pressure drop 

calculations and is given by 

S p = 0.0328 (V s .H s /H t ) (26) 

The total pressure drop is the sum of the two i.e. (_P = A.P + ^ 

OPTIMAL DESIGN 

Optimization, as used here, is the determination of the design of 
an unit process such that it achieves desired effluent goals as for as 
possible and at the same time results in minimum capital and operation 

costs during the design life time of the system. The term "design” 

means the sizes of the system components and their operational 
requirements. There are three preliminary steps that must be completed 
prior to the optimization of a treatment process. First, the system 
must be defined in terms of unit process model. Second, specific input 
values of parameters of the model must be determined. Third, 
limitations on operating characteristics of equipment must be 
obtained. With these sets of information and appropriate cost data, 
optimization is simply a matter of formulating the optimization 
problem and finding a suitable solution technique. 

Unit Process Model 

The process model used to describe ammonia stripping tower 
(counter and cross current), has been described in earlier sections. 
To achieve the desired effluent ammonia concentration, preheating of 
water and/or air, particularly at low operating temperatures, may be 
required. Thus two options are considered: 



1* Design of tower without preheating water and/or air, and 
2. Design of tower with preheating water and/or air. 

Process Data Requirements 

Successful application of the opt iraization procedure described 

here requires that appropriate values of the process data are 

available. These include height of transfer unit for heat transfer HT, 

and height of transfer unit for mass transfer H._. Also the type of 

tG 

packing used should be suitably described (by assigning suitable 
values to H^, V^, and V^) and corresponding value of fanning 

friction factor f should be assumed. These values can be obtained from 
pilot plant studies or the literature (Roesler et al . , 1971). 
Formulation of the Optimisation Problem 

Once the process models and input data have been obtained, the 
next step in system optimization is the formulation of the problem 
itself. In general, optimization problems have the following format. 

Minimise F - f (x^ x £ , x 3 , x r ) (27) 

subject to the constraints 

fijCXi* x 2 . x 3 x n ) 4 0; j*l,2, . . . ,m (28) 

The independent variables of the objective function are called 
decision or design variables. They represent quantities that define 
the system. The set of specific values of these variables which 
minimise the objective function is defined to be the optimal design. 
The set of Eqs . 28 make up the constraint set. These equations set 

limits on the design variables. They are derived from limitations on 
equipment, space or other scarce resources, production requirements, 
or any other restrictions on variables that are functions of the 
design variables directly or indirectly. For ammonia stripping towers, 
the objective function is chosen as the sum of capital (converted to 
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annual) and annual operation costa. Capital coat can be converted to 
annual cost through capital recovery factor (CRF) for a given interest 


rate and system life time. 

Annual Coat * Capital Coat . CRF (29) 

where, CRF can be calculated as 

CRF - [ I Cl + l) n 3/C (1 + I) n - 1] ..(30) 


Here, I is the interest rate, and n is the design life of the system 

or its components. Thus the general form of the objective function is 

TC ■ (Total Capital Coat) CRF + Total Annual Operation Coat 

(31) 

where, TC is the total annual cost. The constraint equations are 
derived from considerations about practical and process aspects. The 
choice of design variables is the first step in formulating the 
objective function. A logical way to proceed is to first determine the 
variables that the cost equations are functions of, and than select 
design variables that fix values of these. 

Counterflow Tower Without Preheating Water and/or Air 

In this case the overall objective function can be written as 

TC « C l + C 2 + C 3 + C 4 + C 5 + C 6 + C 7 + C 8 + C 9 C32) 

where, C^ is the cost of fan, is the cost of pumps, C^ is the cost 

of the tower structure, C. is the exterior and electrical cost, C r is 

4 5 

the cost of the distribution system, is the cost of packing 

material, is the fan power cost, Cg is the pump power cost, and 

is the cost of chemical, labour and supplies. 

An examination of the process model equations presented in 
earlier sections reveal that the designer has the choice of only two , 
variables, namely G' and L’ . Once these values are fixed, all other 
process design variable are uniquely determined. Thus the opimization 



20 


of counterflow tower involves a design vector with two variables G* 
and L ’ . Constraints on these variables were obtained from literature 
(EPA, 1975) and practical considerations. Mathematically, 


1000 4 L* « 20000 

1.0 4 G'/L’ 4 8.0 

Z >, 3.0 


• C33) 


(34) 

(35) 


The lower limit on liquid loading rate is put to ensure that all 
portions of the tower are equally wetted. The higher limit is put so 
that no flooding occurs in the tower. 

Count erf low Tower with Preheating of Water and/or Air 

In this case the overall objective function can be written as 


1 c i ,c i u i ,c t ,c i ,c i ,c i ,c «* c ) 

+ C 10 * c u 


C3<) 


where, C^ is the cost of heating air, and C^ is the cost of heating 
water. As mentioned earlier, the efficiency of ammonia stripping 
process reduces drastically at lower temperature. This warrants 
investigating the option of preheating water and/or air. Thus the cost 
of heating air and water must be incorporated into the cost of 
stripping tower. Therefore the optimization in this case involves a 


design vector with four variables G ' , L ’ , the temperature of air after 


heating Ta^ , and temperature of water after heating Tw^ . The 


additional constriants on Ta^ and Tw^ can be expressed as 


Tw x l Ta x (37) 

Tv 1 l Tw (38) 

Ta^ >, Ta (39) 

Tw x « 70.0 (40) 
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Crossflow Tovar Without Preheating Water and/or Air 

In this case the overall objective function is same as that given 
by Eq. 32. The design vector consists of three variables G', L’ and U. 
The configuration of the tower considered is shown in Fig. 6. The 
constraint in addition to those given by Eqs. 33 to 35 is as follows. 

B/W 4 4.0 ..(41) 

Crossflow Tower with Preheating of Water and/or Air 

In this case the overall objective function is same as that given 
by Eq. 36. The design vector consists of five variables G’ , L ’ , U, 
Ta^ and Tw^. The constraints are given by Eqs. 33 to 41. 

Cost Data 

Cost data is one of the most important input to any optimization 
problem. The cost functions for the ammonia stripping towers were 
taken from Roesler et al (1971). A summary of these cost functions is 
presented in Table 2. The cost for heating water and air for various 
temperature ranges are given in Tables 3 and 4 respectively. These 
values were obtained by considering that water and air were heated in 
heat exchanger with steam as the heating fluid. The cost of heating 
consists of cost of steam, capital cost of the heat exchanger and cost 
of pumping both process and utility fluids through the heat exchanger 
(Peters and Timmerhaus, 1979). In this regard it should be noted that 
the main expense incurred during preheating air and/or water is that 
of the utility fluid i.e. steam in this case. It is assumed that steam 
is available cheaply (as is the case in most industries), and the cost 
of steam was assumed accordingly. It is generally accepted that unless 
sufficient quantity of utility fluid is available cheaply, preheating 
of air and/or water is not feasible. 



Table 2. Summary of cost data for counterflow and crossflow ammonia stripping towers, 

(Adapted from Hoesler et al,, 1971) 
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Table 3. 

Cost in 

cents of 

heating 

1000 Kg of water. 


Initial 
Jemp. , 

10 

Final Temperature , 

20 30 40 

°C 

50 

60 

70 

0 

3.88 

7.44 

10.84 

14.14 

17.38 

20.57 

23.72 

10 

0.0 

3.57 

6.98 

10.30 

13.55 

16.75 

19.92 

20 

— 

0.0 

3.91 

6.74 

9.99 

13.20 

16.38 

30 

— 

— 

0.0 

3.32 

6.58 

9.80 

12.98 

40 

— 

— 

— 

0.0 

3.26 

6.48 

9.66 

50 

— 

— 


— 

0.0 

3.22 

6.40 

60 

— 


— 

— 

— 

0.0 

3.18 


Table 4 

. Cost in 

cents of 

heating 

1000 Kg of 

air. 


Initial 

T emp . , 

°C 

10 

20 

Final 

30 

Temperature, °C 

40 50 

60 

70 

0 

1.85 

3.36 

4.69 

5.91 

7.05 

8.14 

9.19 

10 

0.0 

1.51 

2.84 

4.06 

5.20 

6.29 

7.34 

20 


0.0 

1.33 

2.55 

3.69 

4.78 

5.82 

30 

— 

— 

0.0 

1.22 

2.36 

3.49 

4.49 

40 

— 


— 

0.0 

1.14 

2.23 

3.28 

50 

— 

— 

— 

— 

0.0 

1.09 

2.13 

60 

— 

— 

— 

— 


0.0 

1.05 

The 

capital 

costs of 

structures , fans 

and pumps 

were trans f ormed 


to an annual basis by using Eq. 31 taking design period to be 20 years 
and the rate of interest of 5 percent. This gives a value of CRF equal 
to 0.08024. Further the costs of the ammonia stripping tower and those 
of heating water and air were for different base year and thus they 
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were brought to the same year, 1987, by using suitable cost index 
factors (Engineering News Record, 1987). 

Optimisation Algorithm 

Since the cost functions and constraints are nonlinear in design 
variables , interior penalty function method originally proposed by 
Caroll (1964) has been used. In this method the constrained nonlinear 
optimization problem is converted to unconstrained problem. A new 
function ( - function) is constructed by augmenting a penalty term 

to the objective function. The penalty term is chosen such that its 
value will be small at points away from the constraint boundaries and 
will tend to infinity as the constraint boundaries are approached. 
Caroll (1964) proposed the penalty function as 

9 . 

( X, r ) = F( X ) - r T [ 1/g ( X ) ] (42) 

j = l 3 

where, r is the penalty parameter. The penalty terra in £q. 42 is not 

defined if X is infeasible i.e. any of the g ^ (X) > 0. It has been 

shown (Fiacco and McCormic , 1968) that as r — > 0 

rain (fit X , r ) — > F ( X*) 

* * 

where, F (X ) is the minimum of F (X) and X is the design vector at 
the optimum point. The penalty function is sequentially minimized for 
decreasing sequence of r, until r approaches zero. This is known as 
Sequential Unconstrained Minimization Technique (SUMT) and uses any 
(e.g. Davidon-Fletcher-Powell method with cubic interpolation for one 
dimensional minimization as used here) of the unconstrained 
minimization methods (Rao, 1978). 

RESULTS AND DISCUSSIONS 

Application of optimization procedures developed in the present 
work is demonstrated through evaluation of optimal design for several 




combinations of influent an# process variables for both counter and 
cross flow packed bed towers with and without heating influent water 
and/or air- The typical results are presented in graphical /tabular 
form in this section. 

Fig. 10 shows a typical plot of objective function and penalty 
function versus the number of SUHT cycles for counterflow tower. The 



Fig. 10. Variation of objective and penalty function with SUMT 
cycles for a counterflow tower . 
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penalty function is always greater than the objective function by the 
positive penalty t era. As the number of cycles increases or penalty 
parameter approaches zero, the penalty term decreases and hence the 
difference between penalty function and objective function narrows 
down. At the end as r > 0, the penalty and objective function are 
essentially the same. It can be seen from Fig. 10 that two different 
starting designs lead to same optimal design indicating independence 
of the optimal result from the starting design vector . Similar results 
were obtained for crossflow towers. 

The design vector for counterflow towers consists of only two 
variables and hence optimization results can be represented 
graphically in two dimensional space (Fig. 11). Feasibl e design space 
is mapped as a subset of the design space (set of all possible 
combinations of design variables, G' and L’) enclosed by constraint 
boundaries (Eqs. 33 to 35). To get a better picture of the nature of 
the objective function surface very close to the optimal point , a 
three dimensional representation of a portion of the feasible design 
space is shown in Fig. 12. 

A detailed study was conducted to investigate the effect of 
influent water and air temperatures and flow rate on the total unit 
cost of counter and cross flow towers . A summary of results is 
presented in Tables 5 and 6. Fig. 13 shows typical variation of cost 
and height of the counterflow tower as a function of influent air and 
water temperatures. The variation in cost components as a function of 
flow rate for crossflow towers is shown in Fig. 14. In general 
following inferences can be drawn from the results presented in Figs. 
13 and 14, and Tables 5 and 6 . 



, 1000 kg/hr/m 
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Fig. 11. Graphical represe r ation of constraints and objective 
function for counterflow towers. 0 



Fig. 12. Representation of the nature of the objective functio 




TibU 5» Summary of optiMl dtiiffti of counttrflo* tovors for difftront iopot condition*. 
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^9 

- The optimal unit cost of the crossflow tower is lower than that 
counterflow towers in most cases under similar input, output 
system parameters. 



of 

and 


- The optimal cost of the packed bed towers (cross and counter 
current) increases with decrease in influent water and air 
temperatures. The rate of increase of cost also increases with 
decrease in temperatures. 


- The major portion of the cost consists of (i) the cost of electrical 
power for pumping liquid and blowing air , and ( i i ) chemicals and 
labour . 


Investigations were also conducted for studying the effect of 
influent ammonia concentrations on the cost . It was noted that the 
change in initial ammonia concentration did not alter the cost of the 
unit significantly. Studies were conducted to ascertain the 
feasibility of preheating water and/or air to improve ammonia removal 
efficiency. The objectives of these studies were to f ind answers to 
the following questions . 

1. Is it possible to achieve any reduction in cost of ammonia removal 
by preheating air and/or water in comparison to towers without heating 
arrangements ? 

2. Is it possible to get the effluent ammonia concentration to below 
prescribed standards so that the effluent could be discharged into the 
receiving bodies without further treatment ? 

Optimal designs were evaluated at different flow rates for 
several sets of parameters such as influent ammonia concentration in 
water, air and water temperatures and a target effluent ammonia 
concentration in water. Fig. 15 shows comparative cost of the units 
at different flow rates, and Figs. 16 and 17 present the various cost 
components i.e. cost of heating water, air, cost of the tower and 
total cost for counter and cross flow towers respectively. It is seen 
that substantial reduction in cost can be achieved by preheating water 




0 2 4 6 8 10 12 

Q, mid 


Fig. 15. Variation of optimal cost with flow 
rate . 



Q, mid 

Fig. 16 Main cost components of the optimum design 
of counterflow tower (with heating). 
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Q, mid 

Fig. 17. Main cost components of the optimum design 
of crossflow tower (with heating). 


and/or air. However, it should be noted that this observation is not 
universal and depends on the cost of heating water and/or air relative 
to that of the tower. 

SUMMARY 

l 

In the present study optimization procedures were developed for 
minimum cost design of counter and cross flow ammonia stripping tower 
with and without the option of preheating water and/or air. The 
procedures are based on mathematical models of packed bed counter and 
cross flow ammonia stripping towers. 

Computer programs were developed for the design and cost estimate 
of both counter and cross flow towers. These programs were then 
suitably modified to include them in the optimizaion program as 
subprograms. Optimization problem was framed as minimization of annual 
capital and operating cost subject to the constraints of process 
requirements. Optimization results of counterflow towers are presented 
in two dimensional space indicating constraint boundaries and cos,t 
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contours in feasible design space. 

Comparison between counter and cross flow towers show that the 
latter are economical in most of the cases studied. The cost of the 
units vary almost linearly with flow rate. It was also noted that the 
cost of the system increases with decrease in influent air and water 
temperatures to achieve a certain degree of efficiency. However, 

influent ammonia concentration does not have any significant effect on 

\ 

the cost of the towers. The question of feasibility of heating water 
and/or air was studied in some detail. It was observed that for a 
given set of input, output and system parameters preheating of water 
and/or air leads to a lower cost. 

Optimization procedures described in this thesis are general, but 
while going through the exercise of problem statement, formulation, 
execution and preparation of manuscript, few points were raised which 
could not be included in the scope of the present work. Following are 
a few suggestions to incorporate some of these points. 

1. The height of transfer unit in a packed bed tower depends on the 
type of packing provided. In the present Investigation wooden slats of 
a particular size were used as packing. It is desired that the user 
should have the option of using wooden slats of other sizes as well as 
other packing materials (e.g. raschig rings). Thus optimal design 
should be worked out for other packing materials. 

2. An attempt should be made to develop an equation to evaluate 
optimal cost from a given set of input, output and system parameters. 
This can be achieved by fitting a regression curve after sufficient 
data have been generated. 

3. The optimal design formulation should also include ammonia recovery 
system using absorption towers. 
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